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Abstract 
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1 Introduction 



This article is part of a research program, the purpose of which is to study nonrelativistic integrable and 
superintegrable systems with spin in real three-dimensional Euclidean space. A recent article pQ was devoted 
to a system of two nonrelativistic particles with spin s = o an< ^ s = 0, respectively. Physically, this can be 
interpreted e.g. as a nucleon-pion interaction or an electron - a particle one. An earlier article [2 was devoted 
to the same problem in two dimensions. 

We recall that a Hamiltonian system with n degrees of freedom is called integrable if it allows n independent 
commuting integrals of motion (including the Hamiltonian) and superintegrable if more then n independent 
integrals exist. 

In this paper, we will consider the Hamiltonian 

H = ~A + V (r) + V 1 (r)a-L, (1.1) 

in real three-dimensional Euclidean space E3. We assume that the scalar potential Vo(r) and the spin-orbital 
one V\(r) depend on the (scalar) distance r only. The same system was already considered in |T] and the search 
for superintegrable systems was restricted to those that allow integrals that are first-order polynomials in the 
momenta. 

Here we will concentrate on the case when the integrals of motion are allowed to be second-order matrix 
polynomials in the momenta. Our notations are the same as in jT , i.e. 

p k = -ihd Xk , L k = -ihe k i m xid Xm , (1.2) 

are the linear and angular momentum respectively and 

vi = ( I I ) , °2 = ( °. ) , a 3 = ( I _ X ) , (1.3) 

are the Pauli matrices. 

In the case of a purely scalar potential 

h =\p 2 + Vo(x), (1.4) 

the same assumption, namely that the potential is spherically symmetric, Vq — V(r) would lead only to two 
superintegrable systems, namely the Kepler-Coulomb one Vq = ^ and the harmonic oscillator, Vq — lot 2 . Indeed, 
Bertrand's theorem [31 0] tells us that the only two spherically symmetric potentials in which all bounded 
trajectories are closed are precisely these two. On the other hand, if a classical Hamiltonian (11.41) is maximally 
superintegrable (2n — 1 functionally independent integrals of motion in E n that are well defined functions on 
phase space) then all bounded trajectories must be closed [S]. 
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The superintegrability of the hydrogen atom in quantum mechanics is due to the existence of the Laplace- 
Runge-Lenz integral of motion. This was implicitly used by Pauli [B] and explicitly by Fock [7] and Bargmann 
[5] to calculate the energy levels and wave functions of the hydrogen atom. Similarly, the superintegrability of 
the isotropic harmonic oscillator is due to the existence of the quadrupole tensor [9 (also known as the Fradkin 
tensor [TU]). 

If the potential Vq(x) in (jl.4l) is not spherically symmetric, many new possibilities occur. The first one 
studied was the anisotropic harmonic oscillator with rational ratio of frequencies [9]. 

A systematic study of superintegrable systems in E n with n — 2,3 and n general was started more then 
40 years ago [H]-[TS]. Most of the earlier work was on second-order superintegrability, i.e. with integrals of 
motion that are second-order polynomials in the momenta. The existence of complete sets of commuting second- 
order integrals of motion is directly related to the separation of variables in the Hamilton- Jacobi or Schrodinger 
equation, respectively |11]-|17). Superintegrability of the system (1 1 .4[) is thus related to multiseparability. 

Second-order superintegrability has been studied in 2- and 3-dimensional spaces of constant and nonconstant 
curvature [H]-[IZ],[I8]-[25] and also in n dimensions [26 , 27J|2S]. For third-order superintegrability see [29|-|34 j . 

Recently, infinite families of classical and quantum systems with integrals of arbitrary order have been 
discovered, shown to be superintegrable and solved |35j-|46j. 

Previous studies of superintegrable systems with spin have been of three types. One describes a particle 
with spin interacting with an external field, e.g. an electromagnetic one |47]-|53j. The second describes a spin ~ 
particle interacting with a dyon [54 or with a self-dual monopoles (55j . The third one is our program to study 
superintegrability in a system of two particles of which at least one has nonzero spin [TJ [2] . 

In this paper, we look for superintegrable systems of the form (jl.ip . The system is integrable by construction. 
Since Vo(r) and Vx(r) depend on the distance r alone and (<?, L) is a scalar, the Hamiltonian H commutes with 
the total angular momentum J. As a matter of fact, trivial integrals (present for arbitrary functions Va(r) and 
V\(r)) are 



In view of the rotational and parity invariance of the Hamiltonian (jl.ip , we shall search for integrals of motion 
that have a well defined behavior under these transformations. Thus, we will separately look for systems that 
allow integrals that are scalars, pseudoscalars, vectors or axial vectors. Tensor integrals are left for a future 
study. 

To obtain nontrivial results, we impose from the beginning that the spin-orbital interaction be present 
(Vi 0). We also recall a result from [TJ, namely for 



H, J = L+-a, («r,L), P. 



(1.5) 




.2 ' 



Vq (r) arbitrary 



(1.6) 
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the Hamiltonian (jTTTJ) allows 2 first-order axial vector integrals of motion. They are J and 

S = -~a+H^(x,a). (1.7) 

For 

*(r) = £, %(r) = £, (1.8) 
it allows 2 first-order axial vector integrals and 1 first-order vector integral. They are J, £ and 

n=p-^(xhff). (i.9) 

These two systems are first-order superintegrable and the term V\ (r) — can be induced from a Hamiltonian 
with V\(r) = by a gauge transformation [T]. 

2 0(3) multiplets of integrals of motion 

Rotations and reflections in E(3) leave the Hamiltonian ([l.ip invariant but can transform the integrals of motion 
into new invariants. Thus instead of solving the whole set of determining equations we simplify the problem by 
classifying the integrals of motion into irreducible 0(3) multiplets. 

At our disposal are two vectors x and p and one pseudovector a. The integrals we are considering can involve 
at most second-order powers of p and first-order powers of a, but arbitrary powers of x. 

We shall construct scalars, pseudo-scalars, vectors and axial vectors in the space 

{{£}" xpx a}. (2.1) 

The quantities x, p and a allow us to define 6 independent "directions" in the direct product of the Euclidean 
space and the spin one, namely 

(x, p, L — x A p, a, a A x, a A p\ , (2-2) 

and any 0(3) tensor can be expressed in terms of these. The positive integer n in (|2.1I) is arbitrary and any 
scalar in x space will be written as f(r) where / is an arbitrary function of r = \/ x 2 + y 2 + z 2 . Since p figures 
at most quadratically and a at most linearly we can form exactly seven linearly independent scalars and seven 
pseudoscalars out of the quantities (|2.2I) : 
Scalars: 

Si = l, S 2 =p 2 , S 3 = (x,p), S4 — (a, L) , 

S 5 = (x,p)(a,L), S 6 = L 2 , S 7 = (x,p) 2 . (2.3) 

Pseudoscalars: 

Pi = {a,p) 1 P 2 = (a,x), P 3 =p 2 (x,a), Pi = (x,p)(x,a), 

P 5 = (x,p) 2 (x,a), P e = (x,p)(p,a), P 7 = (x,a)L 2 . (2.4) 
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The independent vectors and axial vectors are as follows. 
Vectors: 

V\ = x , V 2 — p , V 3 — x A a , V4 — p A 3 , V 5 — p 2 x , 

V 6 = p 2 (x A a) , V 7 = (x,p)x, V 8 = {x,p)p, % = (x ,p) (x A 3) , 

Vio = (x,p) {p A 3) , Vn = (3,L)x, V12 = (3, L)p, V13 = {x,p} 2 x , 

V u = (x,p) 2 {x Act), V 15 = (x,p)(a,L)x, V 16 = L 2 x , V 17 = L 2 (x A 3) , 

V ls = (3,x)L, V 19 = (x,p){x,3)L. (2.5) 

Axial vectors: 

Ai = 3, A 2 = L, A 3 =p 2 a, Ai = {x,p)3, A 5 = (x,p)L, A e = (3,L)L, 

A 7 = (x,p) 2 a , A & = (3,p)S, A 9 = (3.p)p, A 10 = (3,x')x, A 11 = (3,x)p, 

A12 = p 2 (x, a) x , A13 = (x,p) (3,x) x , A 14 = {x,p) 2 (x, 3) x , A 15 = (x, p) (p, 3) x , 

Ai & = L 2 3, A 17 = {x,3)p, A 1& = (x, a)L 2 x. (2.6) 

An arbitrary function f(r) is also a scalar and each of the quantities in (|2.3[) - (|2.6[) can be multiplied by f(r) 
without changing its properties under rotations or reflections. 

Even though all the above expressions are linearly independent, higher order polynomial relations between 
them exist. More importantly, we shall use linear relations with coefficients depending on the distance r that 
exist, namely 

S 7 = r 2 S 2 -S 6 , P 7 = r 2 P 3 -P 5l V 16 = r 2 V 5 -V 13 , V 17 = r 2 V 6 - V u , 
V 18 = -r 2 V 4 + V 9 + V n , V 19 = -r 2 V 10 + V u + V 15 , A 16 =r 2 A 3 - A 7 , 

A 17 = r 2 (A 9 - A 3 ) + A 6 + A 7 + A 12 - A 15 , A 18 = r 2 A 12 - A u ■ (2.7) 

We mention that the vector (3,p)L was eliminated from the list (|2.5[) using the nontrivial linear relation 
(a,p)L = V& — Vio + Vi2- We use the relations (|2.7j) to remove the left hand sides of (|2.7p from the analysis 
completely 

The relations determining S 7 , P 7 , Vie, V17, Aiq and Aig are all consequences of the simple vector relation 
L 2 = x 2 p 2 — (x,p) 2 . Those determining Vis and V19 follow from the identity (3,x)(x A p) = — x 2 (p A (?) + 
(x, p) (x A 3) + (3, x A p)x. 

3 Symmetrization of the integrals of motion 

In the rest of the article we separately take the linear combinations of all the scalars, pseudo-scalars, vectors 
and axial vectors with coefficients fi(r) that are functions of r. However, instead of having the bare linear 



combinations of these tensors, the full symmetric forms are needed for the analysis of the commutation relations. 
Thus, in this section, we briefly describe how this symmetrization process carried out basically by working on 
the linear combination of the scalars chosen as a prototype. 

In writing the most general scalar operator, one needs to symmetrize the linear combination of the scalars 
given in (I2.3P 

6 

*fif = £/i(r)£i, (3-1) 

3=1 

term by term. It is obvious that the terms, for example, XiPi and PiXi are in fact different scalars and their 
symmetric form is h{%iPi + Pi%i)- However, at this stage an immediate question arises: should we associate a 
single arbitrary function of r multiplying each term obtained from symmetrization or different weightings are 
necessary for each of them. 

Let us consider the scalar S3, for example. The most general symmetrized form of it is achieved by giving 
an arbitrary function of r to each permutation 

S3 = f a (r)xiPi + f b (r)piXi + Xif c (r)pi + XiPif d (r) + Pif e {r)xi +PiXif f (r) . (3.2) 

Requiring that the operator be hermitian, that is S3 =5*3, we get 

fa{r)xiPi + f b {r)p i x l + Xif c (r)pi + XiPif d (r) +Pif e (r)xi + PiXiff(r) 
= PiXif a (r) + XiPif b {r) + Pif c (r)xi + f d (r)piXi + Xif e (r)pi + f f (r)xiPi . (3.3) 

Hence self-adjointness reduces the half of the arbitrary functions 

fa(r) = ff(r), h{r) = f d {r), f c (r) = f e (r). (3.4) 

Thus, 5*3 can now be rewritten as 

S3 = fa(r)xiPi + PiXif a (r) + f b (r)p l x l + XiPif b (r) + x l f c (r)p l + p l f c (r)x l , (3.5) 

and, using the commutator [xi, f(r)] = 0, we reduce S3 to 

S3 = (fa(r) + f c (r))xiPi +p t Xi(f a (r) + f c (r)) + f b {r)piXi + XiPif b (r) . (3.6) 

Furthermore it is always possible to move the derivative terms to the right 

S 3 = 2(/„(r) + f b (r) + f c (r))x iPi - 3«(/ a (r) + f b (r) + / c (r)) 

-inr(f' a (r) + fl(r) + f c (r)). (3.7) 

Defining 

Mr) = f a (r) + Mr) + f c (r) (3.8) 



we write 

S 3 = 2f 3 (r)x lPl - mh(r) - ihrf'^r) . (3.9) 

Thus for the scalar S3 it is explicitly shown that a single arbitrary function of r is enough to make the 
symmetric form as general as can be. 

However, it is easy to see how difficult this process would become even for the scalars such as S5, for which we 
would have to start by writing 5! = 120 arbitrary functions of r. We have already shown that a single arbitrary 
function of r is enough to give the symmetric form of the operators that are first-order in the momentum pQ . 

To investigate the operators that are second-order in the momentum, we wrote a computer program code 
working under Mathematica. The code organizes the permutations of a given operator into Hermitian couples, 
which does the equivalent steps given in the previous example (j3.3H3.5p . Hence for an operator made of n terms, 
we are left with ^ Hermitian couples. Moving the derivative terms to the extreme right for each couple, we see 
that all of the Hermitian couples for a given operator have the same form (as in (|3.7p ) although some of them 
give a few extra terms. These extra terms can always be absorbed into already existing terms (of lower order 
in p). Hence a redefinition of the arbitrary functions associated with a given operator (as done in Q3.8p ) assures 
that it is enough to have only one arbitrary function of r for each scalar. 

For example, let us consider the symmetrization of the scalar SV There are two distinct Hermitian couples 
associated with S2 , namely: 

Pif2(r)p i +p i f 2 (r)pi (3.10) 
ViPih{r) + h{r)piPi ■ (3.11) 
Moving the derivatives to the right in the first couple, we obtain: 

Pif2(r) Pi +Pif2(r)pi = 2(/ 2 (r)p? - ihf^r) Pi ) . (3.12) 
The same process with the second couple gives: 

PiPiMr) + h{r)PiPi = 2(/ 2 (r)p? - ihf^r) Pi ) - h 2 f^'(r) . (3.13) 

Here, we have an extra term in the second derivative of /2(f)- However, this is just another arbitrary 
function of r and it is already associated with the scalar S%. Hence we could eliminate it by redefining the 
function associated with Si as; fi(r) — ft(r) — ft jf^'M an( i the sum (|3.ip would remain unchanged. The 
Mathematica computer program code makes it clear that the extra terms coming from choosing one Hermitian 
couple over another for a given operator could always be eliminated in this way. 

When writing the operators using tensor notation, one understands that the number of free indices is irrelevant 
to the symmetrization of a given operator, assuring the validity of the above method for any class of operators 
(pseudoscalar, vector and axial vector). In its final version, the computer program is able to take a list of 
operators and symmetrize all of them in an optimal manner keeping no unnecessary terms. 
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4 Commutativity conditions for scalars and pseudoscalars 

4.1 Scalars 

Let us take a linear combination of the independent scalars given in (12.31) 

6 

*s = £/i(r)Si, (4.1) 

2=1 

and fully symmetrize it as described in Section [3] The symmetric form of (|4.ip is written as 

X s = h - 3ihf 3 - ihrfe + 2/ a $ + 2/ B (s, (a, L) p) + /e (£, £) 

+ (/ 4 - Sz^/s - ifir/^) (3, L) + (2/ 3 - ^fC) (x,p) . (4.2) 

The requirement that [H, Xs] — 0, gives us the determining equations for this case. The determining 
equations, obtained by equating the coefficients of the third-order terms to zero in the commutativity equation, 
become 

/b = 0, /£ = 0, $ = 0. (4.3) 

The determining equations, obtained by equating the coefficients of the second-order and first-order terms to 
zero in the commutativity equation, read respectively 

fa = , hfi = 4/ a V{ , (4.4) 

and 

K 2 f[= 4/2^ '. (4.5) 

The rest of the determining equations are then satisfied identically. It is immediately seen that the only solutions 
for equations (|4.3I) - (|4.5p are 

fi=4cxV + c 3 , h = c x , f 4 =4:C 1 V 1 +c 4 , fa = C2, (4.6) 

where Cj, i = 1, . . . , 4 are real constants. 

Thus the corresponding four integrals of motion are those given in 11.51 i.e. there are no nontrivial scalar 
integrals of motion. 

4.2 Pseudoscalars 

As an integral of motion we take a linear combination of the independent pseudoscalars given in (|2.4I) 

6 
2=1 



8 



which has the following fully symmetric form 

X P = (a, f) (/ a - ifi(±/{ + 4/ 4 + r/i)) + (2/ 1 - ift(2/ 3 + 4/ 6 + r/£) ) (<?, £) 
+ (a, f) (2/4 - + IO/5 + 2r/£ + ^)) (x,p) + 2/ 6 (a?,p) 

+ (CT,f)(2/ 3 (p,p)+2/ 5 (f,(i,p)p)) . (4.8) 

Requiring that the commutator [iJ, Xp] = 0, we obtain the determining equations for this case. Those, 
obtained by equating the coefficients of the third-order terms to zero in the commutativity equation, read 

2r/ 5 Vi + ft/ 5 = , (4.9) 

ft/J + 2r (fi/ 5 + (/ a + / 6 ) V^) = , (4.10) 

(ftr + 4r 3 ^) / B + ft/£ - r (2/ 6 Vi - 3ftr/ 5 ) = , (4.11) 

(4r 2 Vi - ft) rft/ 5 - hft + hf' & ~r (2(/ 3 + 2/ 6 ) Vi - 3ftr/^) - , (4.12) 

ft/3 + ft/e + r (ft/3 + 2r (ft/ 5 + / 6 Vi)) = , (4.13) 

h& + 2hf 6 + r (2(/ 3 - h) Vi + 2/5 (2ft + r 2 V 1 )+ 3ftr/ 5 ) = , (4.14) 

ft/3 + ft/e + ^ (2/3 Vi + / B (3ft + 2r 2 Vi) + 2ftr/^) = , (4.15) 

h + U + r (3r/ 5 + f' z + r 2 f' b + = , (4.16) 

ft/ 6 + / 3 (ft-2r 2 Vi) =0, (4.17) 

(2r 2 Vi - ft) (/ 3 + r 2 / 5 + / 6 ) - rft/^ = . (4.18) 

Equation (|4.16p can immediately be integrated to give 

h = Cl - ri r f * + fe) , (4-19) 
where c\ is a real constant. If we introduce (|4.19[) into f|4. 18[) we obtain 

fa = a ( 2Vi - 4") ■ ( 4 - 2 °) 



Then writing fe from (|4.17p and using the rest of the equations (|4.9[) - (|4.15[) . as a compatibility condition we 
obtain a nonlinear second-order equation for Vi 



3ft 4 r(ft - 2r 2 Vi) V{' + Ah 2 [3h 3 + r 2 Vi [9h 2 + 5r 2 Vi (2r 2 Vi - 3ft) J J V{ + 10ft 4 r 3 V/ 2 

+2rV"! 2 (45ft 4 - 4r 2 Vi (l5ft 3 + 4r 4 Vf (r 2 Vi - 3ft))) = . (4.21) 

Performing a standard symmetry analysis |56j for the equation (|4.2ip . we find that the symmetry algebra is 
spanned by two vector fields 

vx = rd r - 2V 1 d Vl , v 2 = r 3 d r + ( J - ^Vx) d Vl , (4.22) 



with [ui,^] = 2«2- It is now possible to lower the order of the equation (by two) using the standard method 
of symmetry reduction for ordinary differential equations [56 . In the case of equation (|4.21[) . this leads to an 
implicit (general) solution which we do not find useful. 

An alternative is to use the symmetry algebra to find particular solutions of (|4.21l) . invariant under the 
subgroup generated by v\, or that generated by v%. The corresponding subgroup invariants are the potentials 

respectively. Substituting (14.231) into (I4.21[) . we see that Vi is a solution for 

Ci = {-|,0,5a|}. 

The invariant V* is a solution only for C2 = and that solution is already included in (I4.24[) (Ci = §)■ The 
solutions V\ for C\ as in (|4.24l) (with C\ ^ $) can be extended to one-parameter classes of solutions by acting 
on them with the symmetry group generated by V2 ■ 

^=(^-^)d-2Ar»)* + A (1 _ 2 V») 1 |A|<^. (4.25) 

Substituting V\ (r) = % into (|4. 25[) and expressing V\ in terms of f , we obtain 

1 ( 2Ci - ft \ , 
TO = -, ft + „ =? • ( 4 - 26 ) 



2r 2 V Vl + 2Af 2 
The 4 values of Ci ^ § in (|4T24) lead to 4 new potentials 

2e 

2r 2 ' Vl + 2Ar 2 ) ' ' iV ' 2r 2 V~ ' Vl + 2Ar 2 



*M = 53 U + < W = ^(l + ^™), (4-27) 



ft h 3h 
2 ' 2 ' 2 



with e 2 = 1. Thus we have 7 potentials to consider: those in (I4.27[) and the original Vi = with C± ■ 
(Ci = is trivial and C\ = 1 is gauge induced and was considered in pQ). Taking A = in (|4.27l) we recover the 
original cases with C\ = — |, (and Ci = 0, 1), but we prefer to treat them separately. 
Case I: Vi = A 

For this type of potential, ([4. 17ft and (14. 19ft immediately imply 

/ 6 = 0, and / 5 = £|-^, (4.28) 

where c\ is an integration constant. Also the set of determining equations given in (|4.9I) - (|4.18I) together with 
(14.19)) give us 

h =--, (4-29) 
r 

where ci is an integration constant. Then the determining equations, obtained from lower order terms, provide 
us with 

4ciV + c 4 c 3 

/i = -r c 3 , / 2 = , k = — , 4.30) 

r r 
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where C3 and C4 are integration constants. For these values of fj (for j = 1, . . . , 6) all the determining equations, 
obtained from the requirement that the commutator [H,Xp] = 0, are satisfied for any Vq — Vo(r). Since we 
have four arbitrary constants and none of them appear in the Hamiltonian, we have four different integrals of 
motion. Two of them are first-order operators and correspond to the ones that were found in [1] , while the other 
two are second-order operators. They are given as, 

Xj. = (4.31) 
r 

X 2 P = -r(a, p) + -(a, x) (x, p)-- (a, x) , (4.32) 
r r 

r3 (a,x) fl \ 2ih -» 2iH 



Xp = 4^^ \^{p,p) + V j + — (a,x) (x,p) - —(a,p) , (4.33) 
X% = (6ih(x,p)-2(x,(x,p)p) +2r 2 (p 1P fj - 2ihr{a,p). (4.34) 

However, there are some polynomial relations (syzygies) between these integrals of motion which could be 



written as 



Xp = 4XpH, (4.35) 



4 = 21 



P ((X P f ft 2 ) , (4.36) 



so no really new integrals are obtained. 
Case II: V x = - A 

For this type of potential, (|4.20p implies 



k = — +c 2 , (4.37) 
r 



where c\ and ci are integration constants. Then the set of determining equations given in (|4.9I) - (|4.18I) gives 

/i=0, / 3 =--, / 4 = 0, / 5 =0, c 2 =0. (4.38) 
r 

Introducing (|4.37l) and (|4.38l) into the determining equations obtained from lower order terms, we find the 
following equations 

M/ 2 -,/ 2 )-^ = 0, ^- ^+^ =0. (4.39) 

Their solutions are 

/ 2 = -g- + r C3 , % = a=£-, (4.40) 

where C3 is an integration constant. Hence the integral of motion Xp depends on two constants c\ and C3. 
However, Vq also depends on a = Thus, we can choose C3 = 4aci and set c\ = 1. The integral of motion 
for this case can be written as 

Xp = +4ar 2_ 2{p ^ + 1 ((f ^ _ ^ (4 41) 

if 



Since Xp is a pseudoscalar operator it also commutes with the components of the total angular momentum 
Ji (i = 1, 2, 3): [H, Xp] = 0, and [Ji,Xp] — 0, as do all the pseudoscalars obtained below. 
Case III: Vi = M 

For this type of potential, the set of determining equations given in (14.9[) - (|4.18p together with (|4.20p give us 

/i = o, / 3 = -- + ^, A = o, 

r 2 



where c\ and C2 are integration constants. If we introduce these integrals of motion into the set of determining 
equations we get ci = 0. Using the determining equations, obtained from lower order terms, we get 

h _ *« + *c l +vs t (4 43) 

and upon introducing (|4.43l) back into the determining equations we obtain a second-order differential equation 
for V 

r 3 (Vo - r Vq) + 15ft 2 = . (4.44) 

Its solution is 



15ft 



2 



y o = ^rW, (4.45) 



where a is an integration constant. Using (|4.45p we obtain fa from (14.43P 

5ft 2 
2r 



/2=ci (4ar--3 ) . (4.46) 



Since there is only one arbitrary constant for this case (a appears in the Hamiltonian) we only have one second- 
order integral of motion, namely 

Xp= K -^-± l-—+4ar 2 ~ — (x,p)-2(p,p) + -(x,(x,p)jf)\ - - ((x,p) - 2ih) (a,p) . (4.47) 
Case IV: Vi = -A 11+ ) , (3 = 2A 



For this type of potential the set of determining equations given in (|4.9|) - (|4. 18[) together with (I4.20P give us 



r r 



ci yl + /3r 2 ci 



. Cl (2 + r 2 /? + £v /l + ^ 2 )- C2 r(l+ £ Vl + /3r 2 ) 

78 - ^ , (4.48) 

where c\ and C2 are integration constants. Introducing (|4.48l) back into the determining equations we get c\ = 
and then rest of the determining equations give us 



f C3 /TTT^ f 2ft 2 c 2 (l + £ VTT^) 

h = J— — , Vo = Wi . (4.49) 

-e+ VI + /3r 2 
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Finally, for this case, we have two arbitrary constants and hence two integrals of motion. One of them is a 
first-order operator and corresponds to the one already found in PQ and the other is a new second-order operator. 
They are given as 

4 = ~V^WM + ^= = i (@,P)-ih), (4.50) 

P -e + VI + j3r 2 v ' 



2 
^2 

2ifr 



(l + e V/1 + /3T 2 ) (ct, x) (h 2 + 3th (x,p) - (x, (x,p)p)^j 

(2 + 6^1 + fir 2 ) + 2(x,p) (a,p) + 2e^/l + 0r 2 {a, x)(p,p). (4.51) 



CaseV:V 1 = ^^l + ^- r j,/3 S 2A 

For this type of potential the set of determining equations given in (|4.9[) - (|4.18[) together with (I4.20[) give us 



2eci y/l + fir 2 a e , 

k = hc 2 , fz = cir J 8+-c 2 V 1 + pH, 

r r 2 

_ 2 ci (2 + r 2 g + 2e y/l + (3 r 2 ) - c 2 r (2 + e y/l + g r 2 ) 

/5 — „ t j (4.52) 

2 r J 

where ci and c 2 are integration constants. Introducing (|4.52l) back into the determining equations we get c 2 = 
and then rest of the determining equations imply 

fx = Q, U = 0, 

-h 2 Cl (1 + 2 r 2 0) (l + iey/l + /3r 2 + Ar 2 p (l + e y/l + /?r 2 )) + 2r 4 (1+/3 r 2 ) c 3 
h = 2r 3 (l + /3r 2 ) 2 ' 

V = j f 7 + 10r 2 /? + 86 (1 + f5r 2 f 2 ) - " Q . (4.53) 

8r 2 (1 + pr 2 ) 2 \ y ' / 4/3(l + /3r 2 ) 

Finally, for this case, we have only one arbitrary constant (03 = ac\) and hence only one second-order 

operator given as 



2r 4 a 



(1 + fir 2 ) - h 2 {l + 2/3 r 2 ) (l + 46^/1 + /3r 2 + Ar 2 f3 (l + e y/l + [3r 2 )^ 



XP 2r3(l + /3r 2 ) 2 (<?,f) 

2i/i 



— 5 + 3/3r 2 +6ey/l + /3r 2 ) (a,x)(x,p) 



2ih 
r 

2 



l + (3r 2 +4ey/l + /3r 2 (a,p) 



- (1 + /3 r 2 ) (a, x ) (p, p) + -5 (2 + r 2 + 2 e ^l + /?r 2 ) (a, jg) (x, (x, p)p) 

- —y/l + Pr 2 (x,p){a,p). (4.54) 
r 

5 Vector integrals of motion 

Let us take a linear combination of the independent vectors given in (12. 5|) 



i=i 
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and fully symmetrize it as described in Section [3] The symmetric form of (|5.1I) is written as 

X v = x(2f 1 -th(4f 7 + ^+rf^ + (2f 11 -ih(Af 15 + ^+rf[ 5 ))(<j,L) 

+ (2f 7 -ih(^ + ^ + Wf 13 + 2rf[ 3 ))(x,p) 

+ 2 (/ 15 (S,p)(a,L) + f 5 p 2 + / 13 (a?, &P)P)j) 

+ (2 (/ 2 - ih(f B + 2/ 8 + T J±) + A2 (a, L) + A [x,p)))p 

+ (x A (?) (2/3 + ih{f n - 4A - y - rA) + 2/ 6 p 2 

+ 2A4 (af, (x,p)p) + (2/ 9 - zft (/k±M + 2r /{ 4 )) (a? 

- (2U-ih(2f 6 + 4f 10 + f 12 + rf[ ) +2Ao(x,p)) (?Ap). (5.2) 

The requirement [if, Xy] = gives us the determining equations for the vectors. The determining equations, 
obtained by equating the coefficients of the third-order terms to zero in the commutativity equation, are 

A + A = o, A = o, As = o, A 4 -A 5 = o, 

Ao + r 2 A 4 + A = o, /6-A2-0, hfi 5 + 2f 6 Vi = o, (5.3) 

together with the following two differential equations 

(A^i)' + 2r/ 6 y i 2 = 0, 

hfl - 2r 2 (/ 6 ^) ' - 4rVi A = . (5.4) 

It is obvious that A = is a solution of the system given in (|5.4[) . For / 6 /0, (|5.4[) implies a compatibility 
condition for V\ which reads as 

QhrV? - 4r 3 V 1 3 + h 2 V{ = . (5.5) 

The above differential equation for V\ had already been considered in [T]. Depending on the solutions of this 
equation we have several cases. The general solution is given by 

* = ^( 1 + 7t?ip)' (5 ' 6) 

where e 2 = 1. Note that V% = ^ and V\ = are special solutions with (e,/3) = (1,0) and (l,oo), respectively. 
The case V\ = is induced by a gauge transformation had already been considered thoroughly in pQ. 

If we introduce the solutions (|5.6p back into the system (|5.4p . we obtain differential equations for A which 
can be solved. Hence, bearing in mind that A — is also a solution of the system (|5.4j) . we have several cases 
to be considered separately. 



14 



5.1 V\ is a solution of (15.51) 
Subcase I: V x = £s 

2r~ 

For this type of potential (|5.4j) implies f§ — c±r. Thus, introducing this together with the relations given in 
(|5.3p into the determining equations obtained by equating the coefficients of the second-order terms to zero in 
the commutativity equation we get 

h = c 3 -r 2 f T , /4 = c 4 -r 2 / 9 , /V^-TT + TTI: /s = -C6, 

2r 2r z 

f» = j, /u = ^ + /9, /i2 = ^+r 3 /r, c 6 = 2c 4 , (5.7) 

where c; (i = 3, . . . , 7) are integration constants. 

The determining equations obtained by equating the coefficients of the first- and zeroth-order terms to zero 
in the commutativity equation, give us 

a(3h 2 - 4r 3 (2V ' + rV ")) = . (5.8) 

Hence we have two possibilities: Either c 4 = 0, or 3ft 2 — 4r 3 (2V ' + rV ") — 0, which has the following solution 

3ft 2 a ,„ „. 

% =**-7' (5 - 9) 

where a is an integration constant and we have set an irrelevant additive constant equal to zero. 
Ii: Vq is given as in (|5.9p 

Upon introducing Vq back into the determining equations coming from first- and zeroth-order terms we have 
c 4 (2ft 2 - 4ar) + c 7 r h 2 c 3 - 2ft 2 cir - 4c 3 ar 1 

/l = ^2 » h = ^2 ' C 5 = -C 3 . (5.10) 

Then all the determining equations are satisfied. We have 4 arbitrary constants (ci, c 3 , c 4 , cj). Thus, we have 
4 integrals of motion which read 

Xlr = -{(xA a),p)p+ — p -—{ x ,p) + -{oAp) + zft 2 — + — 2 (x A a) , (5.11) 
X 2 = 2xp 2 -2(x,p)p~ ^((xAa),p^+2thp-h{aAp)+ih 2{ -^^ + ^{h 2 -2ar), (5.12) 



= —(x-i(x AS)) , (5.13) 

r 2r v ' 

X* = rL(a,p) + yp- y(?Ap)+4 (iH-frfj) . (5.14) 

The integral (|5.13l) is first order and was obtained in [T]. 
I2: Vo unspecified, c 4 = 

From the determining equations obtained by equating the coefficients of the first- and zeroth-order terms to 
zero in the commutativity equation, we get 

c 5 =-c 3 , /i = ^> c 3 (3ft 2 -4r 3 (2^' + rV ")) =0. (5.15) 
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Thus, we must have C3 = 0, because the other possibility gives us the previous potential for Vq. Upon 
introducing back all information into the determining equations we obtain 

h = I , (5-16) 

and then all the determining equations are satisfied for arbitrary values of Vo(r). We have two integrals of 
motion for this case which are given as in (|5.13[) and f|5 . 14|) . It is worth noting that the commutativity condition 
is satisfied for arbitrary scalar potentials Vb (r) . 
Subcase II: Vi = A f 1 + — 



For this type of potential (|5.4p implies /g = C2yl + ft 7 " 2 - Thus, introducing this together with the relations 
given in (|5.3[) into the determining equations obtained by equating the coefficients of the second-order terms to 
zero in the commutativity equation we get 



/2 = c 3 - r / 7 , fi = c 4 - r /n , / 8 = , / 9 = / : 



11 



, c 4 + ec 4 ^TT^ , er 2 (l + < 9r 2 )(c 2; 9-2e(l + /3r 2 )/ 7 ) 

r 1 + e v/1 + /3r 2 + /3r 2 f 2 + /3r 2 + 2ey/l + (3r 2 ) 

where C3 and C4 are integration constants. 

Upon introducing the relations given in (|5.17p into the remaining determining equations, we have 

/i = 0, (5.18) 

together with C2 = and C4 = 0. However, c 2 = reduces fe to zero which will be analyzed below. 

5.2 V\ unspecified, / 6 = 

For this case let us continue to analyze the determining equations obtained by equating the coefficients of the 
second-order terms to zero in the commutativity equation. Introducing fe — together with the relations given 
in (|5.3p into the determining equations we obtain 

h = c 3 - r 2 f r , U = C4 - r 2 f 9 , f 7 = hi (Vi + rVl) , 

/o = 2 2 ^V|5ft ■ ¥n=^9 + /8(^+rW), (5.19) 

together with the relations 

2r(2c 4 + / 8 ) (2r 2 Ui - 3ft) Vj 2 + 2(7 8 (ft 2 + 6r 2 Vi (r 2 Vi - ft)) - ft 2 c 4 ) V/ + r/ 8 (fi - 2r 2 Vi) 2 V/' = , 
2r(2c 4 + / 8 ) (2r 2 Ui - 3ft) V] 2 - 2(/ 8 (2ft 2 + 6r 2 Vi (r 2 ^ - ft)) + ft 2 c 4 ) V? - r/ 8 (ft - 2r 2 Vi) 2 Vf = , 
(2c 4 + h) (-6hrV 2 + 4r 3 U x 3 - h 2 V() = , / 12 (3U/ + rV(') + f{ 2 (Vi + rV/) = , 
fa (3V/ + rVf) = . (5.20) 

Here C3 and C4 are arbitrary constants of integration. From the system (|5 . 20[) we have 4 subcases. 
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Subcase I: /g = — 2c 4 , /12 = cj , 3V[ + rV{' = , which implies 

V 1 = ^+I3. (5.21) 

This is a new spin-orbital potential. 

When we introduce all the information we have found up to now into the remaining determining equations, 
we get 

c 4 (a 2 - (A(3r 2 {/3r 2 - ft))) + r 2 (2ft 2 /1 + 4 c 4 r Vq) = 0, (5.22) 

c 7 = c 3 , c 4 (4 r + 8V ' - 12 p 2 r - ^J^) = . (5.23) 

From equation (|5.23|) it is immediately seen that 

„ , n a(a + 2ft) 
4 r Vq' + 8Vq - 12 /3 2 r - -i— g '- = , (5.24) 

because the other possibility c 4 = 0, either gives us the known potentials (Vi = \ and V\ — A) or reduces the 

number of the integral of motions. The solution of (|5.24[) is given as 

q(q + 2ft)+4/3 2 r 4 - 8r 7 
K, = , (5.25) 

where 7 is an integration constant. Then, introducing back the solution given in (|5.25l) into the remaining 
determining equations, we obtain 

2e 4 (fta-2r( 7 + ft/3r)) c 3 (4r( 7 + gV) - a 2 ) 
h = ^ » /s= 4 r 2 (a + 2/3r 2 ) ' 7 -°' (5 ' 26) 

For these values of fj, (j = 1, . . . , 15) all the determining equations are satisfied for the potentials (|5.21[) and 

(|5.25[) . We have two arbitrary constants, and hence we have two different integrals of motion. They are given as 

„2 



i=x(2p 2 - ( 2/3 ^ 2 (((aA^),p) +ft)) +2(ih-{x,f?))p-h(aAp) 

-^( 2 ^° ) (^f), (5.27) 



V 2r 2 



(l(2ft + a - 2fir 2 ) + ((<? A p+ | (^^) - is) + f (<? A j3) 

ft /2/3r 2 - a\ 

' ^aAJ). (5.28) 



X 2 y 



4 \ r 2 

Subcase II: / 8 = , /12 = , 6hrV? - 4r 3 V? + ft 2 ^' = . 

In this subcase depending on the solutions of (|5.5p we have two possibilities. 

Hi: Vx = ^ 

Since for this type of potential /g in (|5.19p becomes undefined, we need to analyze this case from the beginning. 
Analysis of the determining equations obtained by equating the coefficients of the second-order terms to zero in 
the commutativity equation gives 

/ 2 = c 3) U = -r 2 f 9 , / 7 = 0, /ii = / 9) / 9 = -, (5.29) 

r 
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where C3 and C4 are integration constants. For this case we already have f% = and fa = 0. From the other 
determining equations, we have 

/i = ^> /a = 0, c 3 = 0. (5.30) 

All the determining equations are satisfied for arbitrary potentials Vo (r) and since we have only one arbitrary 
constant, namely C4, there is only one integral of motion for this case. It is given by the equation (|5 . 1 3[) and is 
first-order in the momenta. 

Il2! * = ^( 1 + 7^ 

Upon introducing all the information we have gathered so far into the determining equations obtained from 
equating the first- and zeroth-order terms to zero in the commutativity equation we find 

h -^\ + + J' /3 = ^ ■ (5 ' 31) 

If we introduce the above values of f\ and fs back into the determining equations we see that we must have 
c 3 = and C4 = 0. Then all the determining equations are satisfied. However, if we introduce all the obtained 
values of (i = 1, . . . , 15) into the integral of motion Xy given in (15.2[) it identically vanishes. Hence, we do 
not have any integral of motion for this case. 
Subcase III: 

In this subcase we have /g = 0, C4 = and /12 = 0. Upon introducing this information into the determining 
equations obtained from the first- and zeroth-order terms in the commutativity equation we have 



h = 0, / 3 = ^, V 1 = V — 2 , (5.32) 



where C is an integration constant. Now if we introduce these back into the determining equations we see that 
we either have C3 = or C = 0. If C3 = we have no integral of motion since C4 is already zero. If C = 
then we have V\ = gr? which has been investigated thoroughly. Hence, we conclude that no new information is 
obtained from this case. 
Subcase IV: 

fa (3V/ + rV{') + fa (V, + rV() = (5.33) 

If we solve the above equation for /12 and introduce back into the determining equations obtained from first- 
and zeroth-order terms in the commutativity equation we have 

Ui (Vi + rV()(3V{ + rV(') fa = . (5.34) 

Thus we either have V\ + rV{ = or /{' 2 = 0. The other choice 3V{ + rV{' = has already been investigated. 
V\ + rV{ = implies Vi = — , where C is an integration constant. Then (|5.33[) becomes 

C^=0. (5.35) 
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Hence, we either have C = but then V\ = or f\ 2 = which has been investigated in the previous case. 
The condition /" 2 = implies 

/12 = ci r + en , (5.36) 

where cio and c\\ are integration constants. Introducing this back into the determining equations we see that 
we have cio = and then (|5.33[) becomes 

cu(3Vl + rV(') =0. (5.37) 

The two cases en = and 3V{ + rV{' = have been previously investigated. 
Hence, we conclude that no new information is obtained from this case. 

6 Axial vector integrals of motion 

Let us take a linear combination of the axial vectors given in (I2.6|) 

*A = X>( r )^> (6.1) 
3=1 

and fully symmetrize it as described in Section [3] The symmetric form of (|6.ip is 
X A = 



(a, x)(2f u - ih(2f l2 + U + + 2/ 9 (<?,Pj)p 

+ (2/ 2 - ih{hh + rfi) + 2f 5 (x,p) + 2f 6 (L,af)L + a(2f x - ih(f n + 3/ 4 + fs + rft) 
+ (2/4 - ih(f 15 - 2/ 6 + 8/7 + + 2rf!r)){x,p) + 2f 7 (x, (x,p)p) + 2f 3 p 2 ) 
+ x* (iff, x)(2f w - ih(5f 13 + + f' xl ) + rf[ 3 ) + 2f 12 p 2 + 2f u (x, (x,p)p)) 



+ (2/ 8 - iH(2f 12 + 5/ib + /e + J/9 + rfi B ))(ff,fi 



+ (2/ 13 - ift(l2/ 14 + J(2/{ 2 + /( B ) + 2r/( 4 ))(a,f)(f,p) + 2f 15 (x,p)(a,p)) . (6.2) 
The determining equations obtained by equating the coefficients of the third-order terms to zero in the 
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commutativity equation are 

h = , (6.3) 

2rf 7 + & + r 2 ]'i = , (6.4) 

2rf 7 + f^ + r 2 f^0, (6.5) 

fi/i2+2(/ 9 + /3)Vi=0, (6.6) 

2f 3 V! + 2r 2 /y^i - (/is + /12 + r 2 /i 4 )(^ - 2rVi) - M/is + f' t ) = , (6.7) 

2r/ 14 Vi + ft/( 4 = , (6.8) 

ft/12 + K3S/i4 + 2(/i2 + /is + / 7 )Vi) = , (6.9) 

ft/g + Kft/12 + ft/is + 2/ 3 V a - - , (6.10) 

ft/12 + 2(/ 9 + / 3 )V a + r(Rr/ w + 2r/ 7 ^ + fi/J) = , (6.11) 

2r(/is + r 2 /i4 + 2/ 7 )Vi - ft/( 5 = . (6.12) 

Equation (|6.4[) can be integrated to give 

h = -r 2 h - ci , (6.13) 
where ci is a real constant. Introducing (|6.13p into equation (|6.5[) and integrating, we get 

/ 7 = / 6 + c 2 , (6.14) 
where c 2 is an integration constant. We also solve (|6.6[) for /12 

ft/12 = 2(ci - h + c 2 r 2 + / 6 r 2 )^ . (6.15) 



With the obtained knowledge in equations (|6.13[) - (|6.15[) we solve algebraically for derivatives of the functions. 
Equations (|6.4j) - (l6.12j) are then satisfied for arbitrary values of V\(r). 

Satisfying the determining equations, obtained by equating the coefficients of the third-order terms to zero in 
the commutativity equation, for arbitrary values of Vi(r) we continue with the determining equations, obtained 
by equating the coefficients of the second-order terms to zero in the commutativity equation. We note that four 
of the second-order determining equations involve only the functions / 4 , /§, /n and /13: 

U = 0, (/ii-/a)Vi=0, h(f 11 +U + f a )-2r 2 V 1 f 11 =0, hf[ 1 + hrf 13 + 2rV 1 (f n +U)=0. (6.16) 

They immediately imply that 

/ 4 = 0, / 8 = 0, /n = 0, /i 3 = 0, (6.17) 
which greatly simplifies the rest of the analysis. 
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The complete set of determining equations is too long to present here. Instead, we present a nonlinear 
equation for V\ which is obtained as a compatibility condition for /g in the analysis. Different solutions of this 
condition will form the different cases to be analyzed separately in which the full set of determining equations 
are fully investigated. The compatibility condition reads as 

ci^ = 0, (6.18) 

t 3 

where 

T x = -20H 2 r A V( Vf - 60n 2 r 3 ^ 4 + 16r 7 V? - hh\V[ 2 - m 4 Vi{V( - rV{') , 

T 2 = 3h 5 Vl +r( - Zh^Vl'ih- r 2 V x ) + h 2 r(VxV((-h 2 63 + 20r 2 Vi(3h- r 2 Vi)) - hh\V[ 2 

+ 4rVf(-30ft 3 + r 2 V r i(45n 2 -4r 2 t/i(6ft-r 2 l/i))))) , 
T 3 = -3h 4 r 2 Vl'(2V 1 (3h 2 - 6hr 2 V! + 4r 4 V?) + h 2 rV{) + 6r (Wi^/ (3?i 4 + 4rVi(-9?i 3 

+ 2r 2 V l {m 2 +r 2 V l {bh-2r 2 V l )))) + 2r 3 V 1 4 ( - I5h 4 + &r 4 V^(9h 2 - 2r 2 V 1 (4h - r 2 V^))) 

+ h 4 rV{ 2 (m 2 -2{)r 2 V 1 {h~r 2 V 1 ))y (6.19) 

The compatibility condition (|6.18p is satished if c\ = or I"i = or T 2 = 0. For I\ = or r 2 = a standard 
symmetry analysis is performed in a similar fashion as for equation (|4.21l) and wc find invariant solutions for Vi . 
The symmetry algebra for Ti = is spanned by two vector fields 

vi = rd r - 2V 1 d Vl , v 2 = r 3 d r - Zr 2 Vid Vl , (6.20) 

with [vi,v 2 ] = 2v2- The invariant spin-orbit potentials are respectively 

Vi = ^j and Vl = % , (6.21) 

where C\ and C 2 are integration constants and only for the following values 

C 1 = {-H,~,0,^,H}, and C 2 = {0}, (6.22) 

these two invariants are also solutions of T\ = 0. Thus these values of the constants C\ and C 2 give us the 
special solutions. The group transformations generated by vi are 

Vi = V 1 (l- 2Xr 2 )i , |A|<-^. (6.23) 



VI - 2Ar2 ' 1V ; ' 1 1 2r2 

Upon introducing the value of V\ in the equation (|6.23p . V\ becomes 



f 2 V 1 + 2Af 2 . 

New special solutions of Ti =0 can be obtained from its known solutions, say V\ = with the constant G\ 
taking values from the set given in (I6.22p . since V\(f) is also a solution if Vi(r) is. Hence, we have the following 
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solutions 



c ^- h =* ^ = -f(ttW)*' (6 - 25) 



r 



Cx = =*> 14 = , (6.27) 

« = ft V,^^) 4 . (6.29) 

Together with the four solutions obtained from Vi = ^ with the constant C% taking values from the set given 
in (|6.22p . we have eight different type of potentials V\ for Ty = 0. 
The symmetry algebra for T2 = is spanned by two vector fields 

«i = rd r - 2V 1 d Vl , v 2 = r 3 d r + (ft - 3r 2 Vx)d Vl , (6.30) 
with [wi,U2] = 2^2 . The invariants are 

Vi = and 7i = -f + -| . (6-31) 
where G\ and C2 are integration constants and only for the following values 

C x = {0,^Ay,2ft}, and C 2 - {0} , (6.32) 

these two invariants are also solutions of T2 = 0. Thus these values of the constants Cx and C2 give us the 

special solutions. The group transformations induced by vi are 

~ 1 — 9\r 2 / \ 1 

Vx= r2 [h+(h-r 2 VxWl-2Xr^, |A| < . (6.33) 



Vl - 2Ar 2 

Upon introducing the value of V\ in the equation f|6 . 33[) , Vx becomes 



I / . , ... / ]_ 

1 + 2Af 2 



* = 3 \*+(*-Ci)\ • (6-34) 



New special solutions of T2 = can be obtained from its known solutions, say V\ = with the constant Cx 
taking values from the set given in (I6.32[) . since Vx(f) is also a solution if Vx(r) is. Hence, we have the following 
solutions 

c '= =* ^4( 1 + vfflp)' < 6 - 35 > 

Cx= h Vi = A f (6.37) 

d= ? =► V 1= *( 2 --=i=V (6.38) 

2 2r 2 V Vl + 2Af 2 / V ' 

Cx= 2h =► = 4 f 1 ~ ; 1 V (6-39) 

r 2 l V1 + 2A? 2 / V ' 
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Together with the two new solutions of the form V\ — p>, we have six different type of potentials V\ for this 
case. 

Finally, if we consider the simultaneous solutions of T\ = and T2 = by substituting V" from Ti = in 
T2 = 0, we obtain 

6ftrVi 2 - 4rV 1 3 + h 2 V{ = 0. (6.40) 
Solving this equation provides us with two new potentials 

Vj = A(l + - 6 ), e 2 = l, (6.41) 

where (5 is a constant of integration. 

Although we have 16 different type of potentials which satisfy the compatibility condition (|6.18[) . we note 
that for the following potentials 



2r 2 ' 1 2r 2 ' 1 2rVl + 2Ar 2 ' 

"' = i( 2 + 7rw)- + ,6 ' 42) 

T3 becomes identically zero. Hence we do not consider these potentials in the rest of the analysis. Indeed if 
we consider them and start the investigation of the determining equations from the very beginning the analysis 
gives c\ = 0. This possibility will be analyzed separately. Thus, excluding the potentials given in (|6.42[) as well 
as the potential V\ = which is a gauge induced one and was completely investigated in [T], we now investigate 
the following five cases separately. 
Case I: ci = 

For ci = the third-order determining equations (|6.4j) - (l6.12j) are satisfied by 

f r = Q, / 9 =0, /i2 = 0, /i4 = 0, /i 5 = 0, (6.43) 

together with equations (|6.13[) and (|6.14l) . Then the determining equations, obtained from lower order terms, 
provide us with 

/i = c 3 , / 2 = c 4) /io = 0, C4 = 2c 3 -y, (6.44) 

where C3 and C4 are integration constants. For these values of fj, all the determining equations, obtained from 
the requirement that the commutator [H, Xa] = 0, are satisfied for any Vq = Vo(r) and V\ — Vi(r). Since we 
have two arbitrary constants and none of them are appeared in the Hamiltonian we have the following integrals 
of motion 

X\ = ihUa,S)p + x{5,fj\ -HL- 2(a,L)L- 2itia(x,p) = {(a,L),J) , (6.45) 
X\ = 4L + 2ha = 4 J . (6.46) 
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They are however both trivial since J and (<J,L) are in the set (11.5[) . From here on, we shall list only 
nontrivial integrals. 
Case II: V x = -£ 

For this type of potential, the set of determining equations given in (|6.4[1 - (|6 . 1 2[) . is satisfied by 

/e = -C2, / 9 = 2ci, /i4 = 0, / 15 = -lj, (6.47) 

together with the equations given in (|6.13[) - (|6.15[) . Here, c\ and C2 are integration constants. If we introduce 
these integrals of motion into the set of all determining equations, we get 

/ 2 = c 3 , /io = 0, /i = j + y + c 4 r 2 , (6.48) 

where C3 and C4 are integration constants. Upon introducing all the information obtained so far back into the 
determining equations we find 

V = ar 2 , (6.49) 

with a = — Hence, the integral of motion Xa depends on three constants c\, C2 and C3. The nontrivial 
integral of motion for this case can be written 

„2 , -2\ 



X A = -{2ar z +p z )a + 2(a,p)p 

2 

r 



+ ^lx((a,x)p 2 + 2ih(a,$-2(x^(a,pfj +ih(a(x,p) - (a,x)p)) . (6.50) 



Case III: V x = $ 

For this type of potential, the set of determining equations given in (l6.4D - (|6.12p . is satisfied by 

h = --^ ~ c 2 > h = -2ci , /l4 = ^J, /15 = , (6.51) 

together with the equations given in (|6.13[) - (|6.15[) . Here, c\ and C2 are integration constants. If we introduce 
these integrals of motion into the set of determining equations, obtained from lower order terms, we get 



h = -— + c 3 , fw = -— + c 4 , /l = T + y^y r ' ( 6 - 52 ) 



where C3 and C4 are integration constants. Upon introducing all the information obtained so far back into the 
determining equations we find 

3h 2 

V = — + ar 2 , (6.53) 

with a = J?-. Hence, again the integral of motion Xa depends on three constants c\, C2 and C3. Two of the 
integrals of motion are trivial and the third one is 

X A = (3p 2 - 2ar 2 + 1 (ih(x,p) - (x, (x,p)p)) )S - ^ (h + (L, S))L - 2 ([a,p) - p 

■(^3ihr 2 (a,p) - {a,x)(3h 2 -2r i a + 2r 2 p 2 + I2ih(x,p)) - 4{x, {x,p)p)^j . (6.54) 



2x 
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Case IV: Vi = , (i = 2A, e 2 = 1 

For this type of potential, the set of determining equations given in f|6 .4[) - (|6 . 1 2[) . is satisfied by 



h = -C2 - ^ (l + y 2 + £^1 + , /9 = -2ec lv / l + /3r 2 , / M = 0, As = liflYll+ gg , (6.55) 

together with the equations given in (|6 . 1 3[) - (16. 15|> . Here, c\ and C2 are integration constants. If we introduce 
these integrals of motion into the set of determining equations, obtained from lower order terms, we get 

3 8(l + ar 2 ) 4e \ 



4 2 4r 2 I (l + /3r 2 ) 2 1 + ,3r 2 



2 



/a =cs- ^ (l + 2ar 2 + ^=J== ) , / 10 = 0, (6.56) 



where C3 and a are integration constants. Upon introducing all the information obtained so far back into the 
determining equations we find 

Vo = 8r2(1 ^ r2)2 (4 + ^ - ,^ 2 + 4, (1 + P r 2 ) *) + ^ , (6.57) 

Hence, again the integral of motion depends on three constants ci, C2 and C3. The only nontrivial integral 
of motion is 



X A = \2it? (j + -0= + 1 ± 3£ y+^ j (g, g) - 4feyrW (*,p)j p 

4 ^ 2 + 2«, 2 + -^ + ft (l + f 2 +e VTW) ( L,^ L 
3ft 3 8H(ft 2 + ar 2 ) 4eft 3 2 ( /3 2 



ft 

+ -^r(h 6 + — - — v " ' T ' ^=^= + 8ift 2 ( 1 + -r 2 I (x,p) 

2r 2 ^ (l + /3r 2 ) 2 1 + /3r 2 ^1+^ \ 2 / 

- 8fi ^1 + ^r 2 + e 1/1 + Pr 2 ^j (x, {x,p)p) + Ahr 2 (l + fir 2 + 2e y/l + /3r 2 ) p 2 ^j 

+ % (ih 2 (l + l-r 2 - - 3e^TT^) (a,p) + 4eh^TT^(x,p)(*,p) 

r 2 \ V 2 ^1 + ^r 2 ' 

- 2efi(o ! ,x)-v/l + ^r 2 p 2 ^j . (6.58) 
Case V: V 1 = 4 fl + --f_Y = 2A, 6 2 = 1 



For this type of potential, the set of determining equations given in (|6.4p - (|6.12p . is satisfied by 



h = -C2 



2 Cl 



(l + ^r 2 + e v/1 + /3r 2 ) , / 9 = -2ec iv /l + /?r 2 , 



/i4 = ^ (l + |r 2 + eVl+^ :2 ") , /i5 = 0, (6.59) 



4ci / /3 
— 7- 1 H — » 

together with the equations given in (|6.13[) - (|6.15p . Here, ci and C2 are integration constants. If we introduce 
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these integrals of motion into the set of determining equations, obtained from lower order terms, we get 

c 2 c 3 / r 2 Ua - 6/3 2 ) + r 4 (3/3 3 + - 6/3 2e \ 
h = t + tt - c i — : " " " 



4 2 ^ 4(l + /3r 2 ) 2 ^/l + f3r 2 J 

2ci / e \ ci / 3 15 + 4r 4 a 3e (l + 2/3r 2 ) \ 

where C3 and a are integration constants. Upon introducing all the information obtained so far back into the 
determining equations we find 

3h 2 (4 + 5/3r 2 + 4e (l + /3r 2 ) * ) 

8r 2 (l + /3r 2 ) 2 2/?(l + /3r 2 )' 

Hence, again the integral of motion Xa depends on three constants ci, C2 and C3. The only nontrivial integral 
of motion is 

'2ih 2 ( 5(3 2 e/3r 2 



V o = " TT7^^7T2 - ^7T-1^ • ( 6 - 61 ) 



- fr 2 \— — 2 — + . — -2 1 + ^r 2 + e y/l + /3r 2 ) (x,(x,p)p) 

\ 4(l + /?r 2 ) 2 y/1 + fir 2 J \ 2 v ^ 

+ r 2 (l + /3r 2 + 2eVl + /3r 2 )p 2 j + ^ f (a, x ) (4 (\ + ^r 2 + e v 7 ! + Pr 2 ^j (x,(x,p)p) 



2r 2 



+ «(3 + f r» + ' ^ + 4 f J' ) p) - W (3 + 2ft») + ' < 6 + 7 f a ' ) (g, p)) . (6.62) 

7 Conclusions 

The main results of this article are summed up in Table 1. We list all potentials Vo(r) and V\(r) that have at 
least one integral of motion in addition to the "trivial" ones, i.e. those that exist for all Vo(r) and V\(r) (see 
equation (|1.5t ). 

The results are ordered according to the value of V\(r) (in column 2). In column 3, Vo(r) denotes an arbitrary 
function. Throughout, a and j3 are real arbitrary constants, and we have e = ±1. The nontrivial pseudoscalar, 
vector and axial vector integrals are listed in columns 4,5 and 6. We have included integrals of order and 1 (in 
the momenta) already found in [TJ. Indeed, the potentials in row 1 are first order superintegrable. 

We mention that the nonlinear ordinary differential equations (|4.21l) and I"i = 0, r 2 = with Ti and 
r 2 defined in (|6.19l) . are compatibility conditions for the existence of pseudoscalar and axial vector integrals, 
respectively. We have obtained the general solutions of these equations in implicit form. Since a potential must 
be explicit in order to be useful, we only used the particular explicit solutions obtained by requiring that they 

2G 



Table 1: Superintegrable potentials and their integrals of motion 
No V\ Vo Pscudoscalars Vectors Axial Vectors 

(Em-dsn 



(EMU 

(E2ZD,(E2S1) 



j3j 
(IQ2|) 

be invariant under a subgroup of the symmetry groups of these equations. Whenever possible, we enlarged the 
class of solutions by acting on the subgroup invariant solution with the entire symmetry group of the auxiliary 
equations. 

We consider the most interesting cases to be those with vector or axial vector integrals, specially those that 
can be viewed as a Coulomb atom or harmonic oscillator with a spin orbital interaction. 

A study of the algebras of integrals of motion and of the corresponding solutions of the Pauli-Schrodinger 
equation is in progress, as is a systematic search for integrals that are two index tensors or pseudotensors. 
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h 



h 

77 



h 

27^ 



V (r) 
f? 

63U,(g32) 



3h z _ a 
Sr 2 r 



3 ~h + ar2 <SH) 

4 §s f^ + ar 2 623) 

5 -4 ar 2 



6 2* 3^ + ar 2 



a , a a{a+2h)+if3r 4 - 
27^ + P 8T 2 



1 + ^=4= &m,@m 



9 ^( 1 + 7if^) ra ^ 

10 — — (15371) 



ii 7^1 + 71=) 62D 
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